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Abstract. We develop a theory for the interactions between magnetic impurities in 
nanoscopic systems. The case of impurities in quantum corrals built on the (111) Cu sur- 
' face is analyzed in detail. For elliptical corrals with one impurity, clear magnetic mirages are 

■^j- . obtained. This leads to an enhancement of the inter-impurity interactions when two impuri- 

ties are placed at special points in the corral. We discuss the enhancement of the conduction 
electron response to the local perturbation in other nanoscopic systems. 
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PACS. 75.75.+a - Magnetic properties of nanostructures. 

PACS. 73.22.-f - Electronic structure of nanoscale materials: clusters,nanoparticles, nan- 
otubes, and nanocrystals. 



During the last decade, the physics of nanoscopic and mesoscopic systems has generated an 
increasing interest in the community due to technological advances that allow for controlled 
experiments and the possible developments of new applications, nanoelectronics being one 
example @. 

Coulomb blockade in small metallic or superconducting islands and in quantum dots, the 
similarities between the physics of quantum dots and Kondo impurities ^] or the observation 
of spectroscopic mirages in quantum corrals with Kondo impurities || make evident that in 
' many aspects interactions and correlations play a central role in the behavior of these small 
systems. There are some ingredients that distinguish the physics of nanoscopic systems from 
that of macroscopic samples. One of them, and the most important in the context of the 
present work, is the quantification of one electron levels due to confinement effects. Here we 
discuss the interaction between magnetic impurities mediated by conduction electrons, known 
as the RKKY interaction pj. We analyze this mechanism in nanoscopic systems and show 
that in quantum corrals there is a huge enhancement of the impurity-impurity interaction due 
to confinement effects. 

Quantum corrals are built by positioning atoms, typically transition metal atoms, along 
a closed line on the clean surfaces of noble metals. In recent experiments Manoharan et. 
al. H have built elliptical corrals with Co atoms on the (111) surface of Cu. The Cu (111) 
surface has a band of surface states, orthogonal to the bulk states, which can be represented 
as a two dimensional electron gas confined at the surface. The Fermi level is placed at 
450 meV above the bottom of the surface state band. The atoms forming the corral act 
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Fig. 1 - Spin correlation functions C z between the impurity and conduction electron spins (arbitrary 
units): a) open surface; b) elliptical quantum corral with the impurity at the right focus. 



as scattering centers which tend to confine surface electrons inside the corral. If the corral 
fence where an impenetrable wall electrons inside the corral would be perfectly confined and 
the energy spectrum would consist of a set of delta functions at the energies of the bound 
states. The characteristic energy separation between these states decreases as the size of 
the corral increases. In real systems, there is some leaking of the wave function: electrons 
can tunnel through the fence and the bound states acquire a finite life time. The energy 
spectrum inside the corral (the local density of states) then consists of resonances, the width 
of which increases with increasing energy. This structure can drastically change the response 
of the conduction electrons to a local perturbation. Experiments with different transition 
metal ipurities inside corrals can be made. Among them, only Co impurities clearly show 
spectroscopic evidence of Kondo effect. Here we analyze the case of Fe or Ni impurities and 
describe them with an s-d Hamiltonian. We show that the magnetic response is much more 
pronounced when the impurities are inside the corral rather than on an open surface and in 
some cases a magnetic mirage is observed. To illustrate these effects we anticipate the results 
in fig. [I] with the response of the two-dimensional (2D) surface band to a magnetic impurity. 
When the impurity, interacting via an exchange coupling with the surface electrons, is at a 
free surface (fig. (la)), the usual RKKY oscillations in the correlations between the impurity 
and conduction electron spins are observed. For the 2D case these oscillations decay as the 
inverse square of the distance to the impurity. When the impurity is placed at the focus of 
an elliptical corral with a principal axis of 180A and eccentricity e = 0.75, the polarization 
of the surface states is much larger than for the previous case (fig. (lb), note that the scale 
in figures (la) and (lb) is the same). Moreover, close to the empty focus of the elliptical 
corral there is a substantial maximum which could be interpreted as the magnetic mirage of 
the impurity. These results indicate than under some circumstances two magnetic impurities 
located at special points of the corral could strongly interact despite of the fact that they 
could be up to 100 A apart. The effect is so strong that the impurity-impurity interaction 
at these large distances could be even larger that the interaction of two impurities placed at 
short distances on an open surface. 

These results where obtained by considering the following Hamiltonian: 



H Hcorral ~t~ ^imp 

where H corra i describes the surface electrons in the presence of the corral and is given by 
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Hcorral = V / dr^(r)(-— — V 2 + Vcorral (r) ) l\) a (r) (2) 

* — ' J 2m* 

(7 

here ipt( r ) creates a conduction electron with spin a at coordinate r, m* — 0.3m e is the 
effective mass of the surface electrons and 14 orra /(r) = J2i V( r ~ Rj) is the potential generated 
by the atoms forming the corral that are placed at coordinates R^. We consider typically 
36 atoms in the corral fence as in the experiments of reference (3). The second term in 
Hamiltonian (1) represents a magnetic impurity at coordinate Ri mp that interacts with the 
conduction electron via an exchange interaction J, and in usual notation it reads: 

the operator Si mp describes the impurity spin, the components of a are the Pauli matrices 
and ip^(r) = (^j(r), ^j(r)). For the sd- model considered here, the exchange coupling constant 
J is positive (J ~ lmeV). 

We study diferent geometries: elliptical corrals with eccentricities e = , 0.5 and 0.75, and 
a square corral. To do so, we first calculate the surface-electron propagator in the presence of 
the corral G° (r, r', to) that can be expressed in terms of the free-electron propagator g a {r, r', oj) 
which in 2D is well known ||; in the evaluation of this quantity we have used a high energy 
cutoff. The results inside the ellipse are not very sensitive to the shape of the scattering 
potential V(r) but to the size and shape of the corral, that is to the coordinates Hi of the 
scattering centers. The results presented in this work were obtained with V^r) = a8(r) where 
<5(r) is the Dirac delta function. The electron propagator is given by ||: 

G° (r, r» = g a (r, r', u) + a ^ g a (r, R u lu)[1 - a^g a (Rj , r', w) (4) 

here the matrix g is defined as = g a (Ri, R^, lo). 
The total density of states defined as 

p°( w )=-iW drImG° a (T,r,uj+iO+) (5) 

^" a J corral 

is shown in fig. (2a) and compared with the eigenstates of an isolated ellipse with the same 
shape and dimensions Jt|. For low energies the resonances obtained for the corral are narrow 
and coincide with the eigenvalues of the isolated system. As the energy increases the reso- 
nances become broader and are shifted. This is because for higher energies the fence becomes 
more transparent, when the electron wave length A is smaller than the distance d between scat- 
tering centers, electrons can easily tunnel through the fence. Figure (2a) also shows the density 
of states for a corral of the same dimensions made with 24 scattering centers, that has much 
broader resonances and at the Fermi energy the structure is almost completely washed out. 
For the elliptical corral used in the rest of the paper with e = 0.75, which reproduces one of the 
experimental setups of Ref. (3), the Fermi level lies at one well defined resonance. As we dis- 
cuss below, this coincidence is important and determines the behavior of the system. Finally, 
the local density of states at the Fermi energy p°(r, sf) = — — J2a ^ m ^a( r ^ r i 00 + ^ + )\^=e F 
is shown in figs. (2b) and (2c). 

Now we consider the case of a magnetic impurity inside the corral and analyze the corre- 
lation between the impurity and conduction electron spins. The canonical RKKY calculation 
is perturbative in the exchange interaction J. In the present case there is no singular scat- 
tering as in the Kondo model (J < 0) and perturbation theory gives the correct result. An 
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Fig. 2 - a) Total density of states for the elliptical corral with e = 0.5 bounded by 36 (full line) and 
24 atoms (dotted line). The vertical lines correspond to the energies of a hard wall corral, b) and c): 
top view of the local densities of states at the Fermi energy (450 meV) for elliptical corrals with 36 
atoms and e = 0.5 and 0.75 respectively 



alternative way of doing the calculation is to solve exactly a simpler interaction of the Ising 
form S? i/j' (R.i m p)<j z ip(Ri mp ). We have shown that in small systems the spin-spin correla- 
tions along the z-direction evaluated with the Ising interaction reproduce quantitatively the 
exact correlations calculated with the full isotropic interaction Q. For the Ising-type inter- 
action, the electron propagator G a {r, r',cj) in the presence of the impurity can be evaluated 
by including and extra spin-dependent scattering center in the matrix g of Eq. (^) or can be 
written in terms of G°(r,r',w) as: 

G ff (ry,w) = G^(ry,w) + G°(r ) R 4 ^,w)E imp (w)G^(R iTni ,,r',c I ;) (6) 

where 

JS Z a 

Si mp (o;) = — — — — — — - (7) 

This expression is more appropriate to interpret the results. The spin-spin correlation 
defined as C z (Ri mp ,r) =< Sf mp tp^ (r)a z ip(r) > can be put into the form 

rfw[p T (r,w) -pi(r,w)] (8) 

-oo 

where the spin dependent local densities of states are p a {r, w) = — ^ lmG a (r, r, oj + i0 + ). If 
a second impurity is placed at coordinate R^ mp the effective impurity-impurity coupling is 
proportional to C z (R imp ,'R.' lmp )- 

Figures |l| and || illustrate the behavior of the correlation function. In Fig. (3a) a cut of 
C z (Rj mp ,r) along the principal axis of the elliptical corral with the impurity at one focus 
is shown and compared with the same quantity in the absence of the corral, in fig (3b) the 
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Fig. 3 - a) Spin correlation function C z along the principal axis with the impurity at the left focus for 
the at resonance (thick line) and off resonance (dotted line) situations as compared to the open surface 
correlations (thin line) . b) At resonance C z with the impurity far from the foci; c) off resonance C z 
with the impurity at the focus 



correlation function for the impurity away from the foci is presented. In this last situation 
again, two well pronounced maxima close to the foci are obtained. At very short distances 
C z (Ri mp ,r) depends on the high energy cutoff, however as |Ri mp — r| increases it rapidly 
converges to a stable value. 

Comparing the results of C z (Ri mp ,r) with the local density of states at the Fermi energy 
(fig (2c)), it is clear that the correlation function reproduces the structure of the unperturbed 
density of states p (r,£p). To understand this rather intuitive result, we consider a simple 
limit in which both, the coupling constant J and the resonance widths 7 inside the corral 
are small compared to the energy separation Ae between resonances. Note that for all states 
around a resonance, the wave function inside the corral is essentially of the same form, differing 
only outside the corral. To lowest order, the impurity shifts the position of the resonances 
by 77^0- ~ JSf mp a\(p v (Ri m p)\ 2 . Inserting this result into Eqs. (||) and (||), after some algebra 
and assuming that the Fermi energy coincides with one of the resonant states (at resonance 
situation) we simply obtain C z (Rj mp ,r) = A\(p„ F (r)| 2 where A is a proportionality constant. 
If the Fermi energy lies between two resonant states (off resonance situation), to lowest order 
the spin-spin correlation is zero. The high response of the conduction electrons obtained for 
r — > Rimp is due to higher order corrections. The occurrence of magnetic mirages of impurities 
depends on whether the Fermi energy lies at, or close to, one resonance or between two of 
them, an effect also observed in spectroscopic mirages of Kondo impurities [|| . In quantum 
corrals, since the Fermi energy is fixed by the bulk material and the position of the resonances 
depends on the corral size, the system can be designed to have the Fermi level in any of the 
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Fig. 4 - Spin correlations C z for a circular corral with the impurity at the center (a) and for a square 
corral (b). 



two situations. To show the behavior of the conduction electrons response when the Fermi 
energy is off resonance, instead of changing the corral size, we have shifted the Fermi energy to 
428 meV which corresponds to a valley between two resonances in the local density of states. 
In this case the spin-spin correlation C z (R,i mp ,r) shown in figs. (3a) and (3c) presents small 
amplitude oscillations and an antiferromagnetic correlation with the conduction electrons at 
the empty focus S. 

As in the case of the spectroscopic mirage of Kondo impurities ]g|-[ll|l , the magnetic mirage 
obtained in elliptical corrals is due to a symmetry of the system that generates a local density 
of states at the Fermi level with two pronounced maxima close to the foci. This structure 
enhances the impurity-impurity interaction if the impurities are placed precisely at, or close 
to, these points. In all nanoscopic systems however, the local density of states at e F shows 
some structure that will tend to reinforce the magnetic interactions between impurities even 
if the system geometry is not so special provided that the amplitude |</?„(Ri mp )| 2 of the states 
close to the Fermi level is not too small. In figs. (4a) and (4b) we present the conduction 
electron response for circular and square corrals. A comparison of these results with that of 
the impurity at an open surface (fig (la)), clearly shows the enhancement of the amplitude 
of the oscillations in the spin- spin correlation due to confinement effects. Similar effects are 
expected in small metallic grains. During the last decades a lot of work has been done in small 
metallic clusters with sizes varying from that of a few atoms to macroscopic grains. Magnetic 
impurities in these systems should behave as in the examples presented above. 

Recent advances in STM techniques that allow for the injection of spin-polarized electrons 
JT^ ] could be used for a direct measurement of the effect in quantum corrals. The spin- 
polarized current is obtained by using ferromagnetic Fe or Co tips. A small magnetic field 
B at low temperatures could be used to orient the impurity spin parallel or antiparallcl to 
the tip magnetization and the difference in the current could be measured at each point. If 
the impurity spin is parallel to the magnetization tip, the total current is proportional to 
ap\ (r, £p) + bp^(r, ep) where a and b are characteristic of the tip and are proportional to its 
majority and minority spin densities of states. For the reverse orientation of the impurity spin 
the current is bp-\{r, ep) + apx(r, ep) and the difference is Aj oc (a — &)[p|(r, e F ) — Pi(r, £_f)L 
Although the last bracket is not the spin-spin correlation it has the same structure. In fig. g 
the function Ap = [p-f (r, e F ) — Pi(r, e F )] is shown for two different cases. 

In summary, we have studied the response of conduction electrons with quasi-confined 
(resonant) states to a local perturbation produced by a magnetic impurity. We have analyzed, 
using realistic parameters, the case of impurities in quantum corrals. The main conclusions 
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Fig. 5 - Ap for the impurity at the focus (a) and far from it (b), for the at resonance situation 
are: 

i) If the Fermi energy is at resonance, or close to a resonant state, the response of the 
conduction electrons to the impurity is strongly enhanced. 

ii) If the exchange interaction is small, compared with the typical energy difference of 
resonances, away from the impurity position the spin-spin correlation function reproduces the 
structure of the local density of states at the Fermi energy. 

in) Depending on the particular structure of the system, the effective interaction between 
two impurities could be larger at large distances than at short or intermediate distances. This 
occurs in elliptical corrals where a magnetic mirage is formed due to the particular structure 
of the local density of states at the Fermi energy. 

iv) A direct measurement of the electron response to the impurity could be made by 
injecting spin polarized electrons with a ferromagnetic STM tip. 

v) Similar effects, produced essentially by the confinement of conduction electrons, are 
expected in other systems like small metallic clusters or short carbon nanotubes. 

* * * 
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